Out of equilibrium statistical ensembles for mesoscopic rings coupled to reservoirs 
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FIG. 1: Two terminal geometry for a mesoscopic Aharonov- 
Bohm ring threaded by a magnetic flux <£>. 
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We derive non equilibrium statistical ensembles for a ballistic Aharonov-Bohm loop connected to 
several electrodes connected to reservoirs with different chemical potentials. A striking consequence 
of these non trivial ensembles is the emergence of quantum zero point fluctuations of the persistent 
current around the loop. Detailed predictions for the low frequency noise power are given. 

PACS numbers: 73.23.-b, 73.23.Ra, 72.70.+m 

Persistent currents in isolated rings have been the 
subject of detailed theoretical 0, ^ and experimental 
H @' H fl investigations. They are indeed a probe of 
the many body energy level dependance in the external 
magnetic flux but they are also sensitive to the quantum 
state of the ring. Persistent currents in isolated meso- 
scopic networks have been discussed theoretically in the 
context of a semi-classical model in the diffusive limit , 
and in the context of the scattering approach for ballistic 
conductors Observation of persistent currents in an 
isolated set of connected mesoscopic diffusive rings has 
been reported in recent experiments |J. 

The physics associated to persistent currents in a ring 
connected to external reservoirs has first been investi- 
gated in [Io| where the connection between the spectrum 
of the isolated ring and its transport properties in a two 
terminal geometry has been established. The resulting 
h/e periodic oscillations of the conductance have been ob- 
served experimentally |pr| . Connecting a circuit (which 
plays the role of the system) to external leads will also 
affect the persistent current because of the electron's tun- 
nelling into the leads. For instance, connection to a single 
lead has been considered in 12 to model dissipation in 
the ring. Here, we investigate deeply non-equilibrium sit- 
uations such as a two terminal experiment with a high 
bias voltage. In the non-equilibrium case, the reduced 
density matrix of the ring is a non trivial one, not re- 
ducible to any equilibrium density matrix. This effect 
is non perturbative since it survives even in the limit of 
vanishing coupling to the reservoirs. 
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In this letter, we illustrate this idea in the simple exam- 
ple of a mesoscopic ballistic non interacting ring coupled 
to several electrodes that are connected to external reser- 
voirs of different chemical potentials. We first determine 
the ring's reduced density matrix and show that it can be 
used to recover transport properties in the spirit of [ jlO| . 
Then, we show that a dramatic signature of the non equi- 
librium character of the ring's reduced density matrix is 
the presence of quantum zero point fluctuations of the 
persistent current around the ring. With progresses of 
low noise measurements and nano- fabrication techniques, 



these fluctuations may become observable in a near fu- 
ture. This effect is basically due to the increase of elec- 
tron states that contribute to the noise. Let us mention 
that increases of shot noise based on the same general 
idea have also been discussed in the context of beam ex- 
periments [|l3| and of diffusive contacts with strong elec- 
tron/electron interactions p4| . We end this letter by dis- 
cussing how these non-equilibrium issues can be traced 
within the framework of Keldysh's formalism. 

Similarly to Ref. ||, we divide a connected system 
into two types of regions: a ballistic "intermediate" re- 
gion that contains an Aharonov-Bohm loop and external 
electrodes connected to the Aharonov-Bohm loop (see 
Fig. 0). DC-transport can be controlled through a volt- 
age difference V applied between the left and right ex- 
ternal electrodes and the flux enclosed by the Aharonov- 
Bohm loop. Quantities of physical interest are the trans- 
port current flowing across the mesoscopic device and the 
persistent current flowing around the loop. 

The external electrodes and the intermediate region 
can be described using a tight-binding Hamiltonian for 
free electrons. The coupling between the left /right elec- 
trodes and the intermediate island is described by the 
tunnel Hamiltonian : 

W = t Q ,QC+C Q + t a , a C+C a + t b< /3c£c/3 + t(3. b C^C b . (1) 

Since we are interested in out of equilibrium effects, we 
focus on the eV ^> ksT regime. In most cases, we shall 
set T — OK. This model is valid provided that the tun- 
nelling amplitudes from the electrodes to the small sys- 
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tern are much larger than the inelastic scattering rates in- 
side the small system. In this case, tunnelling to the elec- 
trodes completely determine the non equilibrium state of 
the system and reservoirs. The most transparent picture 
of this stationary non equilibrium state can be obtained 
within the Laudauer formalism. The idea of the scatter- 
ing approach is to find a basis of stationary one-particle 
states (the "in" states) from which the exact non equilib- 
rium many body state of the global system (mesoscopic 
device and leads) can be constructed. The prescription 
is then to fill the "in" states coming from each lead up 
to its chemical potential. 

Since we are interested in the persistent current, a 
quantity specific to the small system, the natural object 
to consider is the one particule reduced density matrix 
restricted to the mesoscopic device. To our knowledge, 
this object has seldom been considered within the con- 
text of mesoscopic physics whereas it has been discussed 
long ago in the context of black hole physics. In the dis- 
cussion of Unruh and Hawking effects [|l5) , the degrees of 
freedom behind the horizon are the counterparts of those 
attached to the leads in the present case. The ring's 
one particle reduced density operator can be obtained by 
eliminating the one particle wave functions in the leads 
from the exact stationary Schrodinger equation. This 
procedure yields an inhomogeneous Schrodinger equation 
on the small system, with an effective Hamiltonian con- 
taining boundary terms representing hopping back and 
forth to the reservoirs and a source term associated to the 
incident electronic flux. Solutions to this equation pro- 
vides the decomposition of "in" states's wave functions 
in the mesoscopic system's Hilbert space. 

This decomposition has a very simple form in the limit 
of very small coupling to the leads. It shows Lorentzian 
resonances when the energy of the "in" state coincides 
with the one particle eigenvalue £\ of the isolated meso- 
scopic system's Hamiltonian. Near resonance, the de- 
composition of the "in" state coming from the left elec- 
trode at energy e is given by: 



' ta,a (A|t 



( £ -£ A ) + flY 



(2) 



Herer A = s ^(|t a ,a<o|A)| 2 + |t /s ,6<6|A)| a ) denotes the to- 
tal escape rate given by the Fermi golden rule for an elec- 
tron in the small system stationary state |A) expressed in 
terms Fermi velocity in the reservoirs vr and the lattice 
spacing I. From the lead's point of view, these resonances 
appear in the transmission coefficient: 



t(e) 
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In the case where the level spacing in the mesoscopic 
system is large compared to the total escape rate T\ of 
an electron from the state |A) to the electrodes, the ring's 
reduced density matrix is diagonal within the basis of 



one particle energy levels of the ring (incoherent reduced 
density matrix). Populations n(e\) 
at zero temperature by: 



c\ c\ ) are given 



" (£A) = 2~*T arCtan {-WIT ) "^f ai ' Ctan {^JT 

(4) 

where Tl._r denote the escape rates in the left (resp. 
right) electrode given by the Fermi golden rule. In gen- 
eral, these escape rates do depend on the energy level 
A. Note that n(e\) is nothing but the average of con- 
tributions of all electrodes connected to the small sys- 
tem, the ponderation being provided by ratios of Fermi 
golden rule's escape rates. In the limit of very small 
tunnelling amplitudes, this distribution can still remain 
non trivial, being a step function determined by the 
ratios of tunnelling amplitudes. The values of these 
non equilibrium populations are easily understood in a 
classical way according to the sequential tunnelling pic- 
ture. For instance, when /il > I-Ir and for energy levels 
Mi* < £a < Ml, the current (1 — n(e\))Ti, from the left 
lead must be equal to the current n(e\)TR to the right 
lead leading to a non equilibrium population Tl/T con- 
sistent with eq. (^). We note that local non equilibrium 
occupation numbers have been discussed and observed 
experimentally in mesoscopic diffusive wires | fl(i| . 

For a non interacting system, the persistent current is 
the average value of the derivative of the ring's Hamilto- 
nian with respect to the magnetic flux. Using the meso- 
scopic system's reduced density matrix and under the 
hypothesis of negligible variation of escape rates for all 
the energy levels which are partially occupied, a simple 
expression can be given in terms of the persistent current 
I(<f),H a ) of a ring at fixed chemical potentiel fi a : 



(•5) 



The current I(<p,jj, a ) can be computed directly in term 
of the single electron eigenenergies e n (4>) for the isolated 
ring as ^2 n de ^^ 0(/z a — £ n (0)) or equivalently using the 
scattering matrix of the ring [ fl7| . Formula (|5|) directly 
shows the influence of non equilibrium populations given 
by eq. (Q). Let us also notice that, in this non equilibrium 
situation, the persistent current cannot be derived using 
the derivative of the average energy with respect to the 
flux. As a function of the external magnetic flux, the 
persistent current still roughly has a sawtooth shape but 
with discontinuities when one particle energy levels cross 
the reservoir's chemical potentials. 

At fixed reservoir's chemical potentials, the experi- 
mental signal is given by the Fourier transform with re- 
spect to the magnetic flux. Let us recall that for an 
isolated ballistic ID ring, the n-th harmonic I n is given 
in terms of the Fermi velocity vf, the ring's perimeter L 
and the number of electrons in the ring N. Introducing 
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h = ev F /L: 



I n = Lx— (-1) 



iJV 



(6) 



In the non equilibrium situation described here, at zero 
temperature, the harmonics are given by : 



InKfJ-a)] = h x — x ( V cos {2nn X { h i a )) ) (7) 
7rn 1 1 ' 



where = y "^2"^ corresponds to an effective flux 

associated with the chemical potential fi. In the par- 
ticular case of a two terminal geometry, the nth har- 
monic is modulated by the bias voltage V. As expected 
on physical grounds, the modulation only appears in the 
asymmetric case. Assuming that eV = ji^ — fip is small 
compared to ^ffi = (s/hl + v//*r)/2, we have for the 
U-dependant part of n-th Fourier harmonic: 



Tr 4. in r \\ ■ ( neVL 
tan (2Tmx{H)) sin 



(8) 

In principle, this could be experimentally observed un- 
der conditions similar to the ones necessary for the ob- 
servation of the persistent current in a single ballistic 
ring. We assume the effective electron gas temperature 
to be smaller than Tivf/L. Tunnelling amplitudes must 
be within the range kBT/Ti <C ^l/r ^ vf/L and the 
voltage bias should satisfy eV ^> Tivf/L. For 10 /an 
diameter ring, the level spacing is typically of 500 mK. A 
temperature within the mK range ensures that thermal 
effects do not suppress the persistent current. Voltages 
above 1 mV should be sufficient to create non equilibrium 
populations for many one particle energy levels. The to- 
tal dc-conductance is given by G — %- — S L ^J? . 

A more striking signature of the non equilibrium re- 
duced density matrix can be found in the zero tempera- 
ture noise of the persistent current. Since the "in" scat- 
tering states are stationary states of the total Hamilto- 
nian, the two time correlation function of the persistent 
current can be exactly computed using our previous ex- 
pressions for ip\(e, a). Its Fourier transform Si p (u>) turns 
out to have a complicated expression showing quantum 
structures above frequencies of order T. For the ring 
considered above, at a dc-resistance of the order of 10 to 
1000 h/e 2 , this frequency is still above 1 Mhz. Measure- 
ments of the persistent current as carried for instance in 
Q are done at a much lower frequency. Therefore, ex- 
perimentalists only access to the zero frequency limit of 
the noise. Sj p (0) turns out to be related to the ensemble 
fluctuations Cp = (Ip) — (Ip) 2 of the persistent current 
by Si p (0) = 2Cp/T. At low frequency, Si P (ui) shows a 
lorentzian behaviour: 



Si p (w) 



2T 



r 2 



C F 



(9) 



The characteristic function of the probability distribution 
of the persistent current is given by a superposition of 
partition noises (j n ((f>) — de n (<j))/d<fi): 



?(k)=l[(l + n(e n (<f > ))(e lk '" {4 ' ) -i 



(10) 



Only the partially populated energy levels contribute to 
the non-zero cumulants and therefore these non-zero cu- 
mulants constitute a true signature of the non equilib- 
rium state of the ring. On the contrary, the persistent 
current of an isolated ring has zero noise at zero temper- 
ature. This opens the possibility of an experimental test 
although measuring the persistent current noise can be 
quite an experimental challenge. 

The total noise is a sum of contributions associated 
to the partially occupied one particle energy levels. For 
eV 3> hvp/L, a continuum spectrum approximation can 
be used. Since all energy levels between the smallest and 
the largest lead's chemical potential contribute to the 
noise, the variation of the Fermi velocity vf(e) might be 
taken into account for explicit evaluations. The variance 
Cp of the persistent current is finally obtained as: 



> , v F (e) de 

t~ 1 -^n(e )(l-ne — 

v F Tivf/L 



(11) 



In particular, for the two terminal geometry, and assum- 
ing constant escape rates and constant Fermi velocity 
over the relevant energy range, one gets: 
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\eV\ TpT! 



* hvp/L T 2 



(12) 



which shows that the fluctuation of the persistent cur- 
rent becomes larger than its average value as soon as the 
voltage is larger than the energy level separation in the 
ring. The three terminal geometry has a richer structure. 
Denoting by /ii,2,3 the chemical potentials in decreasing 
order and y a — T a /T, we get: 



Cf 



el 2 



hvF / L 



{ViViiyi - y 3 ) + V 3 y 3 (l - y 3 ) - V 1 y 1 (l - y x )) 

(13) 



Note that sensitivity with respect to the intermediate 
voltage is enhanced with the asymmetry in the couplings 
to the leads. 

As expected, the noise is also a periodic function of 
the magnetic flux. Its Fourier transform with respect 
to the magnetic flux can easily be obtained using the 
free electron dispersion relation. Using v 2 F {e n {4'))/v F = 
E n ((j))/eF, the zero temperature noise can be related to 
the persistent current around a loop connected to a sin- 
gle lead. Ordering the chemical potentials in increasing 
order, we get: 



f>, f^a)-Ip((t>, Ma+l)) n a ,a+l0- — n a ,a+l) 

(14) 
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where n a . a +i denotes the non equilibrium population 
of all energy levels between chemical potentials fi a and 
fi a +i- In the end, the n-th harmonic of the zero temper- 
ature variance of the persistent current is given by: 

I 2 

C P ( 

(15) 

Let us now connect the simple physical picture just devel- 
oped to some important issues raised within the Keldysh 
approach. The Keldysh formalism is a natural frame- 
work for dealing with non equilibrium physics ]T^ |. It 
provides a way to do systematic perturbation theory in 
the tunnelling amplitudes. The starting point of this per- 
turbation theory is an "initial" density operator which is 
partly fixed by imposing the temperatures and chemical 
potentials of the external reservoirs. But this does not 
by itself determine the initial density operator for the 
mesoscopic ring. 

In order to clarify this point, let us discuss the dc- 
current through system. The current flowing through 
the left contact can be expressed as ]Tg| : 

I a , a = (G< 9 > a - G> a g< a ). (16) 

where the g aa s denote Keldysh's Green functions for the 
isolated left lead whereas the G a . a s are for the system 
with the right lead connected. Using Dyson's equation, 
one gets: G< a = (1 + G R t)g < {\ + tG A ). Naively, equa- 
tion (|l6|)'s r.h.s. could depend on the system's initial 
density operator through Keldysh's Green functions for 
the isolated mesoscopic system. Such terms would intro- 
duce infrared divergences in Keldysh's perturbative ex- 
pansion and of course be present in computations involv- 
ing the persistent current around the ring. But the key 
point is that Dyson's equation for the retarded and ad- 
vanced Green's function provide a way to cure this prob- 
lem in the non interacting case since one can show that 
1 + G R t = G^Xg^) -1 vanishes on the system's eigenen- 
ergies. Those terms can also be cured order by order in 
Keldysh's perturbative expansion by imposing that the 
system's "initial" Green's function g > ' < are given by the 
non trivial one particle reduced density operator com- 
puted in this letter. 

Similar issues have been recently discussed in the con- 
text of quantum dots. For example, in Ref. ^0|, the 
authors ask how to represent the effect of a large volt- 
age bias between the two leads connected to a dot in 
the Kondo regime. This particular problem raises the 
important question of finding the correct starting point 
for Keldysh's perturbative expansion in the coupling to 
the reservoirs. In a recent work pl| , O. Parcollet and 
C. Hooley have pointed out the importance of finding 
the right reduced density operator for the impurity spin 
even in the limit of a small coupling to the leads in order 



to obtain a well defined and correct perturbation expan- 
sion. 

In conclusion, we have shown on a specific example 
that the physics of mesoscopic devices connected to ex- 
ternal leads must be described using non equilibrium sta- 
tistical ensembles even in the limit of vanishing coupling 
to the reservoirs. The example of a mesoscopic ring con- 
nected to two leads at different voltages could provide a 
direct experimental test of these ideas. We have shown 
that non zero cumulants of the persistent current are a 
signature of the non-equilibrium reduced density matrix 
of the mesoscopic ring. These ideas can be extended to 
more complicated circuits using the scattering approach 
on general graphs ^2|. Another issue is to derive the 
non equilibrium state of an interacting system such as a 
Luttinger liquid coupled to leads. 
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